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Abstract. The field theoretic renormalization group and operator product 
expansion are applied to the Kazantsev-Kraichnan kinematic model for the 
magnetohydrodynamic turbulence. The anomalous scaling emerges as a consequence of 
the existence of certain composite fields ("operators") with negative dimensions. The 
anomalous exponents for the correlation functions of arbitrary order are calculated in 
the two- loop approximation (second order of the renormalization-group expansion), 
including the anisotropic sectors. The anomalous scaling and the hierarchy of 
anisotropic contributions become stronger due to those second-order contributions. 

Much attention has been attracted to the problem of intermittency and anomalous 
scaling in developed magnetohydrodynamic (MHD) turbulence; see e.g. P-[9] and 
references therein. It has long been realized that in the so-called Alfvenic regime, 
the MHD turbulence demonstrates the behavior, analogous to that of the ordinary 
fully developed fluid turbulence: cascades of energy from the energy- containing range 
towards smaller scales, where the dissipation effects dominate the dynamics, and self- 
similar (scaling) behavior of the energy spectra in the intermediate (inertial) range. 
However, the intermittent character of the fluctuations in the MHD turbulence is much 
stronger pronounced than that in ordinary fluids. 

The solar wind provides an ideal wind tunnel in which various models and 
approaches to MHD turbulence can be tested [3]-[8]. In solar flares, the highly energetic 
and anisotropic large-scale events coexist with small-scale coherent (singular) structures, 
finally responsible for the dissipation. Thus modeling the way in which the energy is 
transferred along the spectra and then dissipated is a difficult task. The intermittency 
modifies the scaling behavior of the higher-order spectra, leading to anomalous scaling, 
characterized by infinite families of independent exponents. 

A simplified description of that situation was proposed in [2]: the large-scale 
field = UiB^ dominates the dynamics in the distinguished direction n, while the 
activity in the perpendicular plane can be approximated as nearly two-dimensional. 
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This picture allows for reliable numerical simulations which show that the turbulent 
fluctuations tend to organize in rare coherent structures separated by narrow current 
sheets. On the other hand, the observations and simulations show that the scaling 
behavior in the solar wind is more similar to the anomalous scaling of fully developed 
hydrodynamic turbulence, rather than to simple Iroshnikov-Kraichnan scaling suggested 
by two-dimensional picture with the inverse energy cascade; see e.g. the discussion in 
[3]. Thus further analysis of more realistic three-dimensional models is desirable. 

In this paper, we study the inertial-range scaling behavior within the framework 
of a simplified three-dimensional model, known as the Kazantsev-Kraichnan kinematic 
model [in], in which the magnetic field is passive (no feedback on the velocity), while 
the velocity field is modeled by a Gaussian ensemble with prescribed statistics; see also 

mm- 

In spite of their relative simplicity, the models of passive scalar fields advected 
by such "synthetic" velocity ensembles proved to be very interesting because of the 
insight they offer into the origin of intermittency and anomalous scaling in the real fluid 
turbulence on the whole: they reproduce many of the anomalous features of genuine 
turbulent mass or heat transport observed in experiments; see the review paper [17] and 
the literature cited therein. 

Owing to the presence of a new stretching term in the dynamic equation, the 
behavior of the passive vector field appears much richer than that of the scalar field: 
"...there is considerably more life in the large-scale transport of vector quantities" (p. 
232 of Ref. [IB]). Indeed, passive vector fields reveal anomalous scaling already on the 
level of the pair correlation function [HI [I2] . They also develop interesting large-scale 
instabilities that can be interpreted as manifestation of the dynamo effect [IHl dH [IE] . 

In the presence of a constant background field -Bf, the dynamic equation for the 
fluctuating part 9i = 6i{t,x.) of the full magnetic field Bi = B^{ni + 6i) has the form 

dtOi + dk{vkOi - OkVi) = Kd^9i + UkdkVi. (1) 

Here Vi = Vi{t,x.) is the velocity field, k, = c'^/Aira is magnetic diffusion coefficient, 
c is the speed of light, a is the conductivity, is the Laplace operator; summation 
over repeated tensor indices is understood. Equation ([I]) follows from the simplest form 
of Ohm's law for a moving medium j = cr(E + vxB/c) and the Maxwell equations 
neglecting the displacement current; see e.g. [I9]. The last term in the right hand side 
of ([I]) maintains the steady state of the system and is a source of the anisotropy; in 
principle, it can be replaced by an artificial Gaussian noise with appropriate statistics. 

In the real problem, the field Vi satisfies the Navier-Stokes equation with the 
additional Lorentz force term ~ (B x curlB). In the Kazantsev-Kraichnan model 
it obeys a Gaussian distribution with zero mean and correlation function (we consider 
below the incompressible fluid) 

{v,ix)v,ix')) = Do 5{t - t') j^^ ^ p^^.(k) exp[ik ■ (x - x')], (2) 
where -Pii(k) = 5ij — kikj/k"^ is the transverse projector (needed to ensure the 
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incompressibility condition diVi = 0), k = Dq > is an amplitude factor, d is 
the dimensionality of the x space. The parameter ^ can be viewed as a kind of Holder 
exponent, which measures "roughness" of the velocity field [TTj. The IR regularization 
is provided by the cutoff in the integral ([2]) from below at k = m, where m = 1/L is 
the reciprocal of the integral scale L. The anomalous exponents are independent on 
the precise form of the IR regularization; the sharp cutoff is the most convenient choice 
from the calculational viewpoints. The natural interval for the exponent ^ is < < 2, 
when the so-called effective eddy diffusivity has a finite limit for m — t- 0. However, for 
the magnetic field a steady state exists only if ^ < < 2. In the following, we consider 
the physical case d = 3, when = 1 [TT] . 

The model ([1]), ([2]) corresponds to the so-called kinetic regime, in which the effects 
of the magnetic field on the velocity statistics are neglected. In this connection, it is 
worth noting that in the full-scale model of the MHD turbulence such a regime is indeed 
realized in the so-called kinetic fixed point of the renormalization group (RG) equations 
[20] . Various generalizations of this model (finite correlation time, compressibility, more 
general form of the nonlinear terms) were also studied [21] |24]. 

The RG approach to the Kazantsev-Kraichnan model is described in [H] in detail; 
here we only recall the main points. The original stochastic problem ([1]), ([2]) is 
reformulated as a certain field theoretic model. The ultraviolet divergences in the 
corresponding Feynman diagrams manifest themselves as poles in ^. They are removed 
by multiplicative renormalization; as a byproduct of that procedure, differential RG 
equations are derived for various correlation functions. They have an IR attractive 
fixed point. This means that in the IR range Ar 3> 1, where A is the reciprocal of the 
inner (dissipation) length, the correlation functions acquire scaling forms with certain 
critical dimensions of all the fields and parameters F of the model. 

The most important part is played by the critical dimensions A„ ^ associated with 
the irreducible tensor composite fields ( "local composite operators" in the field theoretic 
terminology) built solely of the fields ^ at a single space-time point x = (t, x). They 
have the forms 

F^^i = e,,{x)---e,,{x) {e,{x)e^{x)f + . . . , (3) 

where / < n is the number of the free vector indices and n = / + 2p is the total number 
of the fields 6 entering into the operator; the tensor indices and the argument x of the 
symbol F„ ^ are omitted. The ellipsis stands for the appropriate subtractions involving 
the Kronecker delta symbols, which ensure that the resulting expressions are traceless 
with respect to contraction of any given pair of indices, for example, 6i6j — Sij{6k6k/d) 
and so on. The numbers n and / are even or odd simultaneously. 

The quantities of interest are, in particular, the equal-time pair correlation functions 
of the operators ([3]). For these, solving the corresponding RG equations gives the 
following asymptotic expression 

(F„,Kt,x)Ffc,,(t,xO) ^ (/tA^)"^"^'^/' (^^)-A„,-A,., ^^^^^^_^.(^^) (4) 
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with certain scaling functions ({mr). To simplify the notation, here and below in similar 
expressions we omit the tensor indices and the labels of the functions ({mr). 

The last expression in (jlj) is valid for Ar ^ 1 and arbitrary values of mr. The 
inertial-convective range corresponds to the additional condition that mr <^ 1. The 
forms of the functions ({mr) are not determined by the RG equations themselves; their 
behavior for mr — )■ is studied using Wilson's operator product expansion (OPE). 

According to the OPE, the equal-time product Fi{x)F2{x') of two renormalized 
composite operators at x = (x + x')/2 = const and r = x — x' — )• can be represented 
in the form 

Fi(x)F2(a;')^^C^(r)F(t,x), (5) 

F 

where the functions Cp are the Wilson coefficients, regular in m^, and F are, in 
general, all possible renormalized local composite operators allowed by symmetry. 
More precisely, the operators entering the OPE are those which appear in the 
corresponding Taylor expansions, and also all possible operators that admix to them in 
renormalization. If these operators have additional vector indices, they are contracted 
with the corresponding indices of the coefficients Cp- 

It can always be assumed that the expansion in Eq. (|5]) is made in operators with 
definite critical dimensions Ap. The correlation functions (jl]) are obtained by averaging 
equation of the type ([5]) with the weight expS, where S is the De Dominicis-Janssen 
action functional for our stochastic problem; the quantities (F) appear on the right 
hand sides. Their asymptotic behavior for m — is found from the corresponding RG 
equations and has the form (F) oc m^^ . 

From the expansion ([5]) we therefore find the following asymptotic expression for 
the scaling function ({mr) in the representation (jl]) for mr <^ 1: 

({mr) c^'^^Ap {mr)^'' , (6) 

F 

where the coefficients Ap = Ap{mr) are regular in {mr)"^. 

The feature specific of the models of turbulence is the existence of composite 
operators with negative critical dimensions. Such operators are termed "dangerous," 
because their contributions to the OPE diverge at mr — )■ [25] . 

The dimension of the primary field 6{x) is known exactly, Ag = Ai^q = —1 + ^/2; 
the dimensions of the other operators ([3]) are calculated as infinite series in ^: 

A., = A2e + ASe^ + o(a. (7) 

The dimensions A2,o and A2,2 can be found in a closed form [11] by comparison with the 
exact results for the pair correlation function, obtained within the zero-mode approach 
in [II]-[l3]. For general n and /, the first-order term in was derived in [ll] and has 
the form (for d = 3 and up to the notation): 

Ai} = -n{n + 3)/10 + l{l+l)/5. (8) 
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We have calculated the second term in ([7]), which corresponds to the two-loop 
approximation of the RG. The calculation is rather cumbersome and will be discussed 
elsewhere, and here we only present the final result: 
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700 , ^n(n + 3)-2/(/ + l)J. (9) 

In particular, for the two special cases mentioned above this gives 

A2,o = -^-eV3 + 0(a, A2,2 = e/5 + 7^7375 + 0(a, (10) 

in agreement with the exact results of [II]-[I3]. This confirms validity and mutual 
consistency of the zero-mode and the RG+OPE approaches. 

From the leading-order expression ([8]) it follows that the dimensions ([7]) satisfy the 
hierarchy relations: A„ ; > A„ if Z > which are conveniently expressed as inequalities 
dAn i/dl > 0. This fact, first established in [13], has a deep physical meaning: in the 
presence of large-scale anisotropy, the leading contribution in the inertial-range behavior 
mr -C 1 of the correlation functions like (jl]) is given by the isotropic "shell" (/ = 0). 
The corresponding anomalous exponent is the same as for the purely isotropic case. In 
their turn, the anisotropic contributions give only corrections which vanish for mr — 
the faster the higher the degree of anisotropy / is. This effect gives some quantitative 
support for Kolmogorov's hypothesis of the local isotropy restoration and appears rather 
robust, being observed for the real fluid turbulence and the scalar Kraichnan model 



From it follows that the 0(^^) term in the dimension ([7]) satisfles the same 
inequality, dA^^l/dl ~ (2/ + l)(0.015n + 0.028) > 0, and we conclude that the hierarchy 
of anisotropic contributions becomes stronger when the second-order correction is taken 
into account. 

For the most important case of the scalar operator (|3]) with / = from ([9]) one 
obtains 

A^^l ~ -0.009^2 - 0.052n2 - 0.025n, (11) 

which is negative for all n. Thus the anomalous scaling also becomes stronger 
pronounced when the 0(^^) term is included. In this respect, the magnetic model 
differs drastically from its scalar counterpart, where the higher-order corrections to the 
O(^) term are positive, which eventually leads to the disappearance of the anomalous 
scaling for ,^ — )■ 2; see the discussion in [28] . 

To conclude with, we calculated the anomalous exponents in the Kazantsev- 
Kraichnan kinematic dynamo model in the two-loop order (including the anisotropic 
"shells" in the presence of large-scale anisotropy). We found that both the anomalous 
scaling and the hierarchy of anisotropic contributions become stronger due to the second- 
order corrections to the leading terms. It is interesting to see how these results are 
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affected by compressibility, anisotropy and tlie influence of the magnetic field on the 
velocity field dynamics. This work remains for the future and is partly in progress. 

The authors are indebted to L. Ts. Adzhemyan, Michal Hnatich, Juha Honkonen, 
Andrea Mazzino and Paolo Muratore Ginanneschi for discussions. 
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